#ifndef GAUSS_HPP_
#define GAUSS_HPP_

#include "stdafx.h"

/*! \file gauss.hpp
 * \brief declaration of Gauss points as global variables
 */

extern double xi[51][51]; //!< Legendre Gauss points 
extern double wi[51][51]; //!< Legendre Gauss weights 
extern double yi[9][9];	//!< Jacobi Gauss points
extern double oi[9][9];	//!< Jacobi Gauss weights
extern double Pi; //!< Pi

/*! \fn void init_gauss ()
  \brief has to be called once to initialize the global variables of Gauss points
 */
void init_gauss ();


#endif // GAUSS_HPP_ 
